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Abstract. In the pre-recombination era, cosmological density fluctuations can 
naturally generate magnetic fields through Thomson scatterings. In previous studies, 
only the magnetic field generation from the initially- adiabatic fluctuations has been 
considered. Here we investigate the generation of cosmological magnetic fields sourced 
by the primordial non-adiabatic fluctuations based on the cosmological perturbation 
theory, using the tight-coupling approximations between photon and baryon fluids. It 
is found that the magnetic fields from the non-adiabatic fluctuations can arise at the 
first-order expansion of the tight coupling approximation. This result is in contrast to 
the case of adiabatic initial fluctuations, where the magnetic flelds can be generated 
only at the second-order. In a general case where the primordial density perturbations 
contain small non-adiabatic fluctuations on the top of the dominant adiabatic ones, we 
show that the leading source of magnetic flelds is given by the second-order coupling 
of the adiabatic and non-adiabatic fluctuations. We calculate the power spectrum of 
the generated magnetic flelds when the non-adiabatic fluctuations have a blue power 
spectrum, which has been suggested by recent cosmological observations. 
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1. Introduction 

The magnetic fields in galaxies and in clusters of galaxies are observed and the typical 
strength is about l/xGauss. The origin of these fields is not clear yet [H [21 [3l H] . This 
is one of the significant problem in the modern cosmology. It has been considered that 
very weak large-scale magnetic fields generated in the early universe are amplified by 
the dynamo mechanism after galaxy formation. If the large-scale fields exist in the 
early universe, the cosmic microwave background or the large-scale structure formation 
are affected. We have the upper-limit of the large-scale magnetic field , < 10~^Gauss, 
from CMB observations [5], EJ [71 El [9l [lO] . On the one hand, astrophysical activities are 
considered to vanish information of the primordial magnetic fields near astrophysical 
objects but not the fields in the void regions since the objects rarely exist in void. There 
is an attempt to constrain these fields using the delayed or extended 7-ray emissions 
from distant high-energy astrophysical objects, for example, quasars. We have the lower 
limit > 10^^^ Gauss in the void under some assumptions [HI |12l [131 HH [15] . 

Many mechanisms generating primordial magnetic fields are proposed, such as the 
generation from the second-order density perturbation at the pre-recombination era 
[T6l [T71 [T8l [T9l [201 [2T1 [22] . In this era, protons and electrons move interacting with 
photons by the Thomson scattering. As the mass of electrons is much smaller than 
protons, the cross section of the Thomson scattering for electrons is larger and it causes 
the electric fields. The curly components of these electric fields generate the magnetic 
fields at the pre-recombination era. 

The primordial fluctuations are generated in the very early universe, for example, 
at the time of inflation. They evolve into the anisotropics of the cosmic microwave 
background (CMB) and act as the seed of the large-scale structure. In a simplest single- 
field inflation model, only adiabatic fluctuations are produced. However in other models 
such as a multi-field inflation model or a curvaton model, non-adiabatic fluctuations are 
also produced since a new degree of freedom is added in the system p3]. Actually, 
adiabatic fluctuations are known as the main component of the fluctuations from the 
CMB observations. However, the non-adiabatic fluctuations can still exist as a sub- 
dominant component. For example, in refs. [211 [22], it is indicated that the power 
spectrum of non-adiabatic fluctuations may have a very blue spectrum. On the one 
hand, the possibility of generating vorticity by the non-adiabatic fluctuations is also 
studied recently [261 [271 [30] . 

If the non-adiabatic fluctuations exist, magnetic fields are also generated as is 
expected from the fact that magnetic fields follow similar evolution equations as 
vorticities. In previous studies, initially-adiabatic fluctuations are considered for the 
generation of magnetic fields. Even if one starts from the adiabatic initial conditions, the 
non-adiabatic fluctuations such as differences in motion between photons and baryons 
are generated and they lead to the generations of magnetic field^. In this case, however, 

I There is another spcial mechanism of the magnetic field generation from non-adiabatic fluctuations 
relating cosmic strings [28l [29] . 
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it has been shown that the fields can not be generated at the first order of the tight 
couphng approximation (TCA) for the Thomson scattering and we need to consider the 
second order [191 EI] • This means that the amphtude of the generated magnetic fields 
is suppressed by a tight coupling parameter. 

In this paper, we will show that the fields can be generated at the first-order 
expansion of the TCA if the primordial non-adiabatic fiuctuations exist, and show that 
the leading source of the fields is given by the cross term of the primordial adiabatic and 
primordial non-adiabatic fiuctuations. Hence, there is a possibility that the strength of 
the primordial magnetic fields is stronger than that of the fields only from adiabatic 
fiuctuations if the non-adiabatic fiuctuations have a blue spectrum. Therefore it is 
important to examine the magnetic fields generated by the non-adiabatic fiuctuations. 
The aim of this paper is thus to show explicitly that the primordial non-adiabatic 
fiuctuations can generate the magnetic fields at the first-order expansion of the TCA 
and to obtain and evaluate the power spectrum of the generated magnetic fields. 

The plan of our paper is the following. We show that the non-adiabatic fiuctuations 
generate the magnetic fields at the first-order expansion in the TCA in the next section. 
The calculation is almost the same as the process in ref. [21], but we now include the non- 
adiabatic fiuctuations at initial conditions. In section 3, we derive the power spectrum 
of the generated magnetic fields. We present the results in section 4 and summarize in 
section 5. 

2. Generation of the magnetic fields 

2.1. Evolution equation of the magnetic fields 

We need to consider the second-order cosmological perturbation theory because no 
magnetic fields are generated from scalar-type perturbation at the first order. First 
we define the metric used in this paper. We consider the following metric in the Poisson 
gauge under the fiat-Friedmann spacetime: 



where a and drj = dt/a is scale factor and conformal time. The each metric perturbation 
is expanded as 



where Arabic number represents the order of the cosmological perturbation. 

The basic equations are the equations of motion (EoM) of the matter and the 
induction equation of the magnetic fields. In pre- recombination era, photon(7), 
proton(p) and electron(e) suffuse the universe and interact with each other. Since we 
assume the charge neutrality, the number density for charged particles is the same 
n = rip ~ rig. Then the EoMs for each spices are 



c/s2 ^ q2(^^^) ^_ ^ 20) dif + 2xidT]dx' + (1 - 27^) 6ijdx'dx^] , 



(1) 




(2) 




(3) 
(4) 



VvT(p)j'' = enEi + Kf + /tf^. 



Primordial magnetic fields from the non-adiahatic fluctuations 



4 



eriEi + + K 
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(5) 

where Ei and n'f in the right hand side are the electric fields and collision terms 
[201 l32l [33] . The Thomson scattering between photons and charged particles is 



— K 



PI 



rni^ 4 

-^-aTnp^(l-27^) 
mi 6 



(^(7)» 



(6) 

(7) 



where and Il^-'j are the cross section of the Thomson scattering and anisotropic stress 
of photons. The Coulomb scattering between charged particles is 

Kf = -Kf = -e^n2?7c(l - 27^) {v(^p)i - v^e)i) , (8) 

where rjc is the electric resistivity. We also use the induction equation of the magnetic 
fields dMl) 



(a^E*)' = -e'^'^dj [a (1 + 0) Ek] - {avjEk) 



(9) 



where the prime represents the derivative with conformal time rj. 

Subtracting two EoMs for protons and electrons leads to "the Ohm's law". In 
our case, there are two velocity differences, between photon and baryon 6v^b = — Vb 
and between proton and electron 6vpe = Vp — Ve where Vb is the baryon's velocity of 
the center of mass. We assume that proton and electron are strongly coupled because 
coulomb interaction is strong. It results that the velocity difference between protons 
and electrons is very small compared with one between photons and baryons. Even if 
the electric current does not exist, the magnetic field can be generated. As a result, we 
obtain "the Ohm's law" 

1 - Aar 



E,, 



1 + 3e 



■ap^{l - 27^)A^;, 



(7b)i ) 



where we define 

:= Sv(^b)i 



4 



V^b)j^. 



J 

7 i' 



(10) 



:iii 



P = rrie/mp as the ratio between the mass of proton and electron and neglect 5f(pe)j- If 
the mass of proton and electron is same, /3 = 1, the electric fields are not generated. 
Eq. ffTOl) corresponds to the equation added anisotropic stresses to eq.(27) in ref . [2T] . 

Substituting eq. ffTOj) into the induction equation eq.Q, we obtain the evolution 
equation of the magnetic fields 

1 - Aaj 



(3 3e 



X 



dj {p^Av(^yb)k) + P-ydj{4> - 27^)A^;(^fe)fc + — [p^Vja'^Av(^yb)ky 



(12) 



The equation means that the sources generating the magnetic field are the velocity 
difference between photons and baryons and anisotropic stress of photons. 
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2.2. Analytic formula of the velocity difference Af(^b) 

We use the tight couphng approximation (TCA) for the Thomson scattering so as 
to represent Af (^6)^ by the well-known perturbative quantities. When the dynamical 
timescale is longer than the timescale of the Thomson scattering, we can expand the 
physical quantities under the TCA. At zeroth order in the TCA, all fluid components 
move together with the same velocity fj. We set the photon frame as the "background" 
in the TCA. In other words, physical quantity for zeroth-order TCA is the same as that 
of photon. The velocity for baryons is expanded as 

Vib)i = Vi + v[ll + v[l^^ + ■ ■ ■ , (13) 

where Vi is the common velocity of photons and baryons in the tight coupling limit and 
also the photon's velocity. Since various physical quantities such as Sv(^^h)i, are 
expanded 

n^. = <i' + <? + ■••, (15) 

the quantity At>(^b)i in the eq. f|T2l) is expressed 

Av^,,y = Av[!^ly + ^^1% + ■■■■ (16) 

The equation which determines the velocity difference between baryon and photon 
is provided from the subtraction of the EoM for photon from baryon 

— (^(7)i + Xi) - — {v(h)i + Xi) + {5v{ih)i)' + 4:'H6v^^b)i 



- (0 + 27^) 



— f (7)i + {5v{^b)i) + Al-L5v(^i,)i 



- 5n'Sv^^b)i + + (^(7)^^(7)) - YT^^' {^(p)K) + /^^(-)^^(e)) 

= -ail - 27^)Ay(^fe)i, (17) 
where a is defined as 

a:=^— + (18) 

1 + p Snip 

with R := (3mp(l + /3)n)/(4p^). From the next section, we solve At>(^fe)j under the above 
expansion of the tight coupling approximation. 

2.2.1. Background of the TCA Firstly we note the background in the TCA and define 
the non-adiabatic fluctuations. In the zeroth-order TCA and the first-order cosmological 
perturbation, the energy conservations for photon and baryon, 

^'^'^^ -4(7^«)'- W^), f^y = 3(7^W)'-9.^;^W, (19) 
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imply that 



6jl 



4 
3 



n 



(0) 



(20) 



where the bar represents the quantity in the zeroth-order TCA. In the case of primordial 
non-adiabatic fluctuations, there is a difference between the fluctuations of photons and 
baryons as an initial condition which means 



n 



(0) 



4pf 



+ Cix). 



(21) 



In the above equation, C{x) corresponds to the non-adiabatic fluctuation, which is 
usually called baryon isocurvature mode. This indicates that each fluctuations move 
with the common velocity keeping the difference in density contrast which is caused 
by the non-adiabatic fluctuations. In the adiabatic initial condition, this term becomes 
zero, so photons and baryons move together, as one fluid. As we shall show later C{x) 
plays a central role in generating magnetic fields. 

2.2.2. First-order cosmological perturbation The first-order cosmological perturbation 
is the most simple calculation. In this order, eq. f|T7|) presents 



(n(0))' 



-(0) 

p\ 



n 



(0) 



4 



(22) 



and so we obtain 



(Ai) 
(76)1 



a 



(0) 



(1) 



1 di5p. 



7 



4 



where we use background equations and define 

;i + /32)(l + i?(o))4aaTP 



(0) 

'7 



1 + 



(23) 



(24) 



In this order, the non-adiabatic fluctuations do not contribute because the density 
fluctuations of baryon do not appear. 



2.2.3. Second-order cosmological perturbation The effect of the non-adiabatic 
fluctuations and anisotropic stress appear firstly in this order. Equation f lT7|) becomes 



n 



(0) 



n 



(0) ^(0) 



- (0(1) + 27^(l)) 



(^(0))' 



n 



(0) 
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(0) 



(2) , (2) 

VI ' + Xi 



6pj di6p-y di6p- 



(2)' 
'7 



(25) 



where 



a 



a 



(0) 



1 + i?(o) 
4 + 3i?(o) 



-(0) 



n 



(0) 



(1) 

'7 



4(l + i?(o)) 



+ 



1 + 



a 



a 



(0) 



ad 



1 + 



C'(^). (26) 



Baryon density perturbation, Sn^^\ is in both sides of eq. fl2H]) . however, the non-adiabatic 
term in the left-hand side vanishes because the terms in the curly brackets become 
{5n^^y n^^'^)' and non-adiabatic fluctuations are independent of time at this order from 
eq. (l2T|l . This can be solved for ^V(^^h)i S^^^ 



a 



(0) 



a 



a 



(0) 



,(1) 



ad 



di6p. 



(2) 

'7 



+ 



(7^«)' 



(0) 



ad / 



,(1) 



2 Pt^ J 



The last term is the contribution of the non-adiabatic fluctuations and the anisotropic 
stress. 



(27) 



2.3. Source term of the magnetic fields 



Here we show the generation of the magnetic fields from the non-adiabatic fluctuations. 
We substitute eqs. (l23l) and ( 1271) into eq.(fT2l) and obtain the evolution equation of the 
magnetic fields 



1 - Aar 
1 + 3e 



Pi 



<''-2K('))A.;;;i + ^(.^pf.„fA.!;;i)' 



1 - (5^ Aar 
1 + 3e 



where we use 



a'-pfe'^^ 



2a^H 



(2^ 



a 



(Av 



(Ai) V 

{■yb)k) 



+ nAv 



1 + 



(7b) fc' 



(29) 
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with photon's vorticity w*^^)* which is defined in eq. OA.Sp . Equation fl28p tells us that the 
non-adiabatic fiuctuations contribute as the source at the first-order TCA. This results 
from the fact that the direction of the velocity difference between photon and baryon is 
different from the gradient of the non-adiabatic fluctuations. 

The non-adiabatic fluctuations generate also the vorticity in eq. (|28l) . We take curl 
of a total momentum conservation at flrst-order TCA and obtain 

^ ^ 1 + 2(1 + i?(o))2 1 mk \ ) 

If the non-adiabatic fluctuations do not exist, the vorticity decays with time and the 
magnetic flelds can not be generated at this order. That agree with refs.[l9l [21] 

3. Power spectrum of the generated magnetic fields 

3.1. Fourier transformation of the source term 

We will derive the power spectrum of the magnetic fields generated by the non-adiabatic 
fiuctuations. The source term in eq. (|25]) is 

Since we consider first-order scalar perturbations, the velocity is given by a gradient of 
a scalar function, v. Defining 

vf\x\n) := d,v^'\ 5«(x\r/):=^ (32) 

and 

<5V-(i):=W^)- V), (33) 

then we can write ^v^l^l^^ in Eq.dSI]) as ^v^l^l^. = di5V^^^/a^^\ The Fourier 
transformation of eq. (!3T]) is 

^« = + TTmww^ I "'"'^^ " (34) 

where a hat represents the Fourier transformation. We solve eq.(l30l) for the vorticity to 
obtain 

Therefore the Fourier transformation of the magnetic fields is 



1 + /3 ea3 J a-^ ^ 1 + (3 ea^ Jq 



a~ 



§ Eqs. ([28|) and ((30|) look different from eq.(D7) in ref.[22]. However this comes from the difference of 
the frame. Our frame is the same as the frame in refs.[Tni [H]. There are the details in Appendix D in 
ref.Ea. 
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where p^o* is the present photon energy density and 

= (^ll^i / ^'^(^^ X ^'ClP6vt, , (37) 

' ^ ' '''dv'-^ fdMkxj^^C^'^wi'. (38) 



(27r)3/2 aio) 1 + i?(o) ' i + R{o) 
3.2. Evolution of the fluctuations 

In order to evaluate the source terms, we need to know the evolution of 5V^. In the linear 
order, 5V^ is given by the sum of the adiabatic and non-adiabatic modes. Though both 
the fluctuations should be responsible for the generation of magnetic fields, we, hereafter, 
consider only the adiabatic mode in calculating SV^:. This is because the non-adiabatic 
fiuctuations are expected to be smaller than the adiabatic fiuctuation from observations. 
Accordingly, the equations flM|) and f l35|) tell that the source of the magnetic field is the 
cross coupling of the adiabatic fiuctuation (^Vg) and the non-adiabatic fiuctuation (C*^). 

We consider the radiation-dominated era when the Thomson scattering is strong 
and the tight coupling approximation is good enough. Then we can neglect matter 
perturbations and anisotropic stress in linearized Einstein equations, namely = 

5*-^^ + 5*'^V(pmVpi°'') — and 11^ = 0. The linearized Einstein equations become 

7,K m,k ^ III/ y^), I 

m (iiif)' + nnf) + enf = -AnCa^pfsf , (39) 

{nfy + H^f = -47rGa2pf (1 + ^)4^) , (40) 

- 4'^ = ' (41) 
and the evolution equation of the radiation is 

(5«)' + 3n {cl - w) 6f = il + w) [evf^ + SiTZfy) . (42) 

Here the sound velocity, Cg, and equation of state, w, are given by = w = 1/3. Metric 
perturbations are represented by the one variable tp^}^ := TZf = from (gl]). We 
substitute this into eqs. (l39!) . (HO!) and fH2|) and derive the differential equation in terms 

(4'^)" + ^(4'^)' + ^^'€^ = 0' (43) 

where we used Friedmann equations and H = l/r/ for radiation-dominated era. We 
define y := kr]/\/3 and solve this equation to obtain 

^iJ(^)(k,r^) = 34'\kf-^, (44) 

y 

where t/'q^'* is the primordial adiabatic fiuctuation and ji{y) is spherical Bessel's function 
of order 1, 

ji (y) = \ (sin y-ycosy). (45) 

y2 
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In solving the evolution of the fluctuations, we used eq.( H2]) . At small scale, however, 
we need to solve the Boltzmann equation with anisotropic stress. Then the fluctuations 
dissipate by the Silk damping. The dissipating fluctuation can not generate the magnetic 
fields. So we consider effectively this influences by introducing the cut-off scale later. 
From eqs.( l39|) . ( l40l) and f l4Tl) . the perturbation 6v\:'^ is 

^vf = nv^, - \h = ^,k^^, = lv'k%. (46) 

Using the above equation, we can write the source, S\ as 

S^= f d'p {px kyC^^^l^l_.f{k,p,v), (47) 



where 



f{k,p,r]):= , Drn^^-mJ ^~P1 ' (^^ 



(27r)3/2 2(l + i?(o))a(o)' ^' y 
with y = \k — p\ri/\/3. Furthermore the source, can be written as 



= I d'p {px kycl}^^^l^g{k,p,v), (49) 



,) ^ ^±,J^ f^^mim (50) 



where 



(27r)3/2a(o)2r/(l + i?(o)) 7 'l + i^W y 
with y = \k — plfi/y/S. 

3.3. Power spectrum 

The power spectrum of the generated magnetic fields Pb is defined as 

(^B,{k, V)B*%K, r/)) := '^Pj,{k)6{k - K), (51) 

where bracket means to take the ensemble average. From eq.f l36l) . the left-hand side in 
the above equation is 

X ((a(fc,r/i)^*'(^,^2)) + (s,{k,ri,)S*\K,ri2))) , (52) 

where we neglect correlation between 5* and Q, for example (SQ^ for simplicity. 
We estimate (Si{k,r]i)S*\K ,r]2) 

'Si{k,r],)S*\K,r]2)] 

d^pi j d^p2{k X Pl)^{K X p2yf{k,pi,r]i)f{K,p2,r]2) x 

X (c^'\p-'i)^l^\k-p-,)C*^'\p-2)^f\K-p-'2) 
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d^Pi j d^P2{k X pi)i{K X p2yf{k,pi,rii)f{K,p2,ri2) x 
X [{C^'\p-[)i'i'\k - p-{)) {c<'\p-2)i'f\K - p-2)] 
+ (c^'^p-Dc^'^p-^)) (i>^'\k - p-[)i,f\K ^ p-^)] 

where we used Wick's theorem in the second equahty. The power spectrum of the 
primordial adiabatic fluctuations and non-adiabatic fluctuations is defined respectively 

as 

(4(^)4* W) ■= 2n^Paik)5ik - K) (54) 

and 

'C{k)C*{K)) := 27,^ PUk)6{k-K). (55) 



We assume no correlation between adiabatic and non-adiabatic mode, {Cipoj = 0. 
Hence eq. (!53|) becomes 

(^S,{k,r],)S*\K,r]2)) = {2Ti^f5{k - K) j d^'pik x x p)' x 

X Pna{v)Pa{\k-fh)f\k,p,r,^)f{k,p,r,2) (56) 

and similarly 

(a(^,r/l)^i«(i?,r72)) ={27i^f5{k-K) j d'p{k x p),{K x pY 

X Pnaip)Pa{\k - p\)gik,p,r]i)g{K,p,r]2). (57) 
Finally we obtain the power spectrum of the magnetic fields as 

^PB{k) = [j^^^ i^^rj d'p\kxp\'PUp)Pa{\k-p\) X 

X / dT]i dr]2a^'^a2'^{g{k,p,r]i)g{k,p,r]2) + f{k,p,r]i)f{k,p,r]2)}, (58) 
Jo Jo 

where 



f{k,p,r]):= , f,(n,. -^n^ \k-p\ , (59) 



g{k,p,r]) := - / dr] - , (60) 



(277)3/2 2(1 + i?(o))a(o)' ^1 y 

(2vr)3/2a(0) 2r^(l + i?(0)) J ^iTR(<^) y 

_(0) ^ (l + /3^)(l + i?W)aaTpf ^ ^g^^ 
1 + /3 ' 

= (3mp(l + ^)n(o))/(4p?'^), and y = \k-p\r]/^/3. 
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In the end of this section, we note about the power spectra of the adiabatic and 
non-adiabatic fluctuations. The power spectrum of the adiabatic mode, Pa{k), can be 
well described as a power law 



P^{k) = pA^(fco) (^^ j , (62) 

where /cq = 0.002Mpc~^ is the pivot scale, ^ 2.5 x 10^^, and ni ~ 1.0 from the 
7- year WMAP data |31j. We assume that the non-adiabatic fluctuations have also a 
power law. The power spectrum of the non-adiabatic fluctuations of baryon is 

Pna 

1^ [-or) U j ■ 

where A is represented as the ratio of the cold dark matter (CDM) isocurvature 
fluctuation to the sum of the adiabatic and the isocurvature fluctuations. The factor of 
comes in order to convert the power spectrum of CDM isocurvature mode 
to that of the baryon isocurvature one C{x) [34]. Recent observations indicate a blue 
spectrum with ^2 = 2 — 4 and ^ < 0.1 — 0.01, depending on the observational data used 
[311 1211 125]. In this paper, we adopt A ~ 0.01 and n2 = 4 as a simple example. Since the 
amplitude of the generated magnetic fields is proportional to the root of the amplitude 
of the non-adiabatic fluctuations, B oc y/Ps oc \/A, it is insentive for the amplitude of 
the non-adiabatic fluctuations. Therefore the amplitude of the non-adiabatic fluctuation 
is less important. We give the power spectrum of the generated magnetic fields for the 
different index of the non-adiabatic fluctuations later. 



4. Results 



4.1. Integral region and Results 

We need to perform numerical integrations of eg. (1581) . As a cut-off scale we consider 
the Silk dumping which is the dissipation of fluctuations by a radiation drag. The 
fluctuations whose scale is less than the Silk scale dissipate and do not create the 
magnetic fields any more. 

We take it into account that the Silk scale depends on time ksuk — 4.3 x 
lO^^a-^/^Mpc"! [201 [35]. Furthermore we neglect the fluctuations with larger scale at 
each time than the horizon because large-scale modes do not contribute to the generation 
of the magnetic fields. The integral region for a wave number k is presented as the gray 
trapezoid region in Fig. [H 

We perform the numerical integration for A ^ 0.1,0.01 and n2 — 4 to calculate 
the amplitudes of the magnetic fields at the radiation-matter equality time and show 
the result in Fig. [2l The circle and plus mark represent the contribution of the slip 
term (SS*) for each A and the square and cross mark represent the contribution of the 
vorticity term (QQ*) for each A. From this figure we find that the vorticity term is 
less important at the generation of the magnetic fields at pre-recombination era. The 
physical magnetic field strength is about 10~^° Gauss at IMpc. The generated magnetic 
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Figure 1. The integral region of ea. ([58)) . rjeq is the confornial time which is the 
radiation-matter equality time. The horizon is depend on the time, such as kn oc 77^^ 
and the Silk scale is also depend on the time, such as ksuk ry"^/^. Then the 
maximum and minimum of the wavenumber and the conformal time are determined 
for an arbitrary wavenumber, fc, of the magnetic fields. 
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Figure 2. Calculated physical amplitudes of magnetic fields for several comoving wave 
numbers (fc). The plus mark and cross mark represent the contributions of the slip 
term and vorticity term respectively. Here we set A ~ 0.1 (circle and square mark), 
A ~ 0.01 (plus and cross mark) and n2 = 4. Two lines are the analytically-estimated 
power spectra (see text). 



fields have a blue spectrum k^-^^ and so the larger strength is expected at the smaller 
scales. 

4.2. Simple estimation for arbitrary spectral index n2 

Finally we show a simple estimation of the k-dependence of the magnetic-field spectra. 
For the slip term, the power spectrum of the magnetic fields is Pb ~ k^p'^'^^^\k — 
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p|~'^(J fT]~'^dr])'^. Since J rj^'^fdrj ~ ri'^\k —p\~^, Pb ~ k^p^ri'^/\k —p\^. Since the main 
contribution to the integration for p is from p ~ \k ~ pi ~ /c, we obtain Pb ~ k'^'^^^rj^. 
We consider the silk scale rj oc k'"^/^ and obtain _B ~ ^/Pb behavior of 

the magnetic fields. For 722 = 4, 5 ~ We take the same process about vorticity 

term and obtain B ~ a/P^ ~ /.n2/2-7/6_ p^j, ^2 = 4, i? ~ /j5/6~o.83_ rpj^g deviation from 
this power law at the large scale is caused by the boundary effect from the horizon. This 
fact indicates that the contribution at the Silk scale is very important in the generation 
of the fields at small scales. 

5. Conclusion 

In this paper, we discussed the generation of magnetic fields by non-adiabatic 
fluctuations at the pre-recombination era. Firstly we showed analytically that non- 
adiabatic fluctuations generate magnetic fields using the tight couphng approximation. 
We found that magnetic fields are generated at the first-order expansion of the tight 
coupling approximation. This result should be compared with the case of initially- 
adiabatic fluctuations, where magnetic fields are created only at the second order. 

Secondly we calculated the power spectra of magnetic fields at the radiation-matter 
equality time considering the non-adiabatic fluctuations with a blue spectrum. For 
this we only consider the second order coupling between primordial non-adiabatic and 
adiabatic fluctuations because this coupling is expected to give larger magnetic fields 
than the auto-coupling of non-adiabatic fluctuations. We found that the fields have a 
blue spectrum B oc /c^-^^ for n2 = 4, and the amphtude is about B ~ 10^^° Gauss at the 
comoving scale of IMpc if the maximum amount of non-adiabatic fluctuations allowed 
from CMB observations is considered. Because the spectrum is bluer than those from 
the adiabatic fluctuations found in refs. [Sni ES] , this fields may dominate the primordial 
magnetic fields at small scales. The amplitude would be enough to be amplified to the 
present fields in the galaxy by dynamo mechanism [36], however it is still insufficient to 
explain the fields in void regions recently claimed by [HI [121 [13 EH US] • Clearly both 
detailed observations and theoretical investigations are necessary to understand how the 
large-scale magnetic fields evolve after their generation in voids. 

Because we consider only the first order of the tight coupling approximation, we 
cannot directly extend our analysis to the time of recombination. For a complete 
analysis, we have to solve the time evolution of At'(^fe) numerically. This will be our 
future work. 
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Appendix A. Four velocities, vorticities and electromagnetic fields 

In this appendix, we write down four velocities, vorticities and electromagnetic fields 
used in our paper. From the metric and normalized condition u^m'* = —1, the four 
vector of the each species is 

= - + - + Xi- 27^^;(/)^j, (A.l) 

where is the 3-dimensional velocity. The energy-momentum tensor of each 

spices(/ = p, e, 7) is represented as 

where = p^/3 and Pp^e = 0. We give the useful relation 

v.T(,)o" = -[p'j + m{pi + pi)] + 3(p/ + p,)7^' 

- {pi + Pi)div\i), (A.3) 
^uTi^i/ = [{pi + Pi) iv{i)i + X'd] ' + 4'H(p/ + Pi) {v(i)i + Xi) 
+ (pi + Pi)di{.'P - 0^) + diPi 
- (0 + 27^) I [{pi + pi) Vi^i-),] ' + m {pi + pi) v^i)i } 

- 5(p7 + pi)n'v(^i)i + {pi + pi)dj {v(^i)iV^^j^ , (A.4) 

The vorticity of photons is defined as ujI^ := —}^e'^''f^U(^y)\VuU(^^)p where e^'^P^ is the 
alternative tensor with e^^^^a = ij ^~g. Since the first-order vector perturbation does 
not exist, the leading term of the vorticity is second-order perturbation. 

'a, + X-f ) + ^"S/^S; + + 2K'")] . (A.5) 

where e^^^ is the 3-dimensional flat alternating tensor with e^^^ = 1. 

The Faraday tensor F^y is decomposed to the electric and magnetic fields 

= F^-'u^,^,, = . (A.6) 
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